Photonic families of non-Gaussian entangled states and entanglement criteria for 

continuous-variable systems 
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We consider two classes of non-Gaussian entangled states generated from the product of number 
states with the action of the beamsplitter or the two-mode squeezer. It is shown that, for many of 
these states, the covariance matrix is compatible with the covariance matrix of separable Gaussian 
states and their separability cannot be verified by the measurements of the first and second moments 
of canonical variables. We identify a couple of continuous- variable entanglement criteria with higher 
order moments to verify these non-Gaussian entanglement. 



I. INTRODUCTION 

Gaussian states, such as coherent states, squeezed 
states, and two-mode squeezed states, are of fundamental 
to describe basic nature of quantum mechanical phenom- 
ena. They play a central role in experiments of quantum 
optics and quantum information science Gaussian 
states are characterized by their covariance matrices with 
the first and second order moments of canonical quadra- 
ture variables. Similarly, Gaussian operations are charac- 
terized by the transformation of the covariance matrices, 
and there have been various experimental demonstrations 
of quantum information processing within the Gaussian 
framework 0, Q . In turn, it has been known that non- 
Gaussian resource is necessary to implement a universal 
quantum computation 0-0 ■ This fact strongly moti- 
vates us to study the property of non-Gaussian quantum 
states. 

There have been many works to inv estig ate the prop- 
erties of non-Gaussian wave functions [Tol ITT| and non- 
Gaussian entangled states @, HHlij. An outstanding 
example of no n- Gaussian quantum states is a single pho- 
ton state. Generation method of the single photon state 
and its non-Gaussian properties have been widely investi- 
gated not only for its primitive quantum nature but also 
for potential application in quantum information science 
[2l| . Quantum measurements to manipulate a few num- 
ber of photons have been experimentally demonstrated 
[23 - |25| and such methods are useful to generate non- 
Gaussian entanglement QJ, [H [H H, S3 . 

Entanglement of Gaussian states has been successfully 
described by their covariance matrices 0, H, [28l - f30j . On 
the other hand, there have been many proposed entan- 
glement criteria based on the measurements of hi ghe r or- 
der moments of canonical quadrature variables |3lM38{ . 
Generally, the criteria with higher order moments are 
thought to be powerful to detect non-Gaussian entan- 
glement, however, the criteria based on the covariance 
matrices are sensitive to some of non-Gaussian entan- 
gled states. Thus, it is not so clear in what condition we 
need higher order moments for entanglement verification. 
Hence, it would be valuable to investigate limitations of 
entanglement criterion with the lower order moments of 
the covariance matrix and show the advantage of higher 



order criteria [16|, [l7[ . Thereby, the following two steps 
would be important: (i) To identify a set of non-Gaussian 
entangled states whose inseparability cannot be verified 
by the measurement of the covariance matrix, (ii) To 
find higher order entanglement criterion which can verify 
the inseparability of these states. 

The photonic number states are eigenstates of a har- 
monic oscillator and their properties are simply described 
by using the annihilation and creation operators. Such 
a formalism enables us to describe an orthonormal ba- 
sis consist of entangled states [20| . It is likely that such 
an orthonormal basis gives an insightful example to pro- 
ceed the step (i). To proceed the step (ii), there have 
been a sequence of works to present the separable condi- 
tions based on SU(2) and SU(1,1) commutation relations 
[32l [33l I3514371 ] . These conditions can be derived from 
the uncertainty relations for the operators described by 
the products of canonical variables, and are thought to 
be the lowest order criteria to go beyond the Gaussian 
framework. 

In this paper, we consider two classes of non-Gaussian 
entangled states generated from the product of the num- 
ber states by using two-mode Gaussian unitary interac- 
tions. We show that, for many of these states, the covari- 
ance matrix is compatible with the covariance matrix of 
separable Gaussian states and their inseparability cannot 
be verified by the measurements of the first and second 
moments of canonical variables. We also identify a couple 
of continuous-variable entanglement criteria with higher 
order moments to verify the inseparability of these states. 

This paper is organized as follows. We investigate the 
property of the covariance matrix for the non-Gaussian 
entangled states in Sec. II. We consider the entanglement 
criteria with higher order moments for the non-Gaussian 
entangled states in Sec. III. The results are summarized 
in Sec IV. 



II. PHOTONIC NON-GAUSSIAN ENTANGLED 
STATES AND ENTANGLEMENT CRITERIA 
FOR GAUSSIAN STATES 

In this section we apply the separable condition based 
on the covariance matrix of the density operator to two 
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families of non-Gaussian entangled states, and show lim- 
itations on the detection of entanglement from the mea- 
surement of the covariance matrix. 



Covariance matrix and inseparability of 
two-mode states 



We consider a two-mode system AB described by 
quadrature variables with the canonical commutation re- 
lations [i/tjj'ii] = [xb,Pb] = i- Let us define a vector 
form for the set of quadrature variables as 



R := (xa,Pa,xb,PbY- 



(1) 



The covariance matrix of a two-mode density operator p 
is defined by 



7: = (RR l + (RR t ) t ) -2(R)(R l ) 
= ((AR)(AR t ) + [(Ai?)(A#)]*) 



where (O) = Tr(Op) denotes the expectation value of the 
observable O and AO := O— (O). The submatrices A, B, 
and C are real symmetric 2-by-2 matrices. Their deter- 
minants and dot 7 are invariant under the local Gaussian 
unitary transformations. In terms of these determinants, 
the separable criterion [29[ is given by 

D := dct 7 + 1 - dct A — dct B + 2 det C > 0. (3) 

This is the necessary and sufficient condition for separa- 
bility of two-mode Gaussian states. For non-Gaussian 
states, this is a necessary condition for separability. 
Hence, we may define a class of purely non- Gaussian en- 
tangled states by the set of entangled states which fulfill 
the condition of Eq. (J3)). An important property of this 
class of states is that their covariance matrices cannot 
be transform into the covariant matrices of the entan- 
gled Gaussian states by local Gaussian operations and 
classical communications. 

In the next two subsections we consider two families 
of entangled states given by acting Gaussian unitary in- 
teractions on the product of number states as in FIG. [U 
We will specify the regime of parameters which indicates 
the purely non-Gaussian entangled states in each family. 
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(b) 
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FIG. 1: Two families of photonic non-Gaussian entangled 
states from a product of number states, (a) The two-mode 
squeezed number states are generated by the action of the 
two-mode squeezing operation, (b) The beamsplitted number 
states are generated by the beam-splitter operation. The pa- 
rameters £ and r specify the action of the two-mode squeezer 
and the action of the beamsplitter, respectively (see main 
text). 



vacuum state in the Schmidt decomposed form is defined 
by 



(4) 



n=0 



where we assume |£| < 1. The TMS vacuum state is a 
Gaussian state and it is known that any two-mode pure 
Gaussian state can be transformed into the form of Eq. 
(j4]) by local Gaussian unitary operations. 

Let us define the annihilation operators crossing over 
the separation of A and B by 



A s := V s aVj = 



a — ± + b — fa^ , . 

= , B ( = VtW£ = ^-i_ ,(5) 



where the unitary operator := e^ atfcl ^ ab is the two- 
mode squeezing operator. The new field operators fulfill 
the bosonic commutation relations [A%, At] = [B^,Bl] = 



1 



and [A^,B 5 ] = = 0. The TMS vacuum state 

corresponds to the vacuum of the new field operators as 
we have if#f> = and B^) = 0. The TMS number 
state [13, H, H(| is defined by the number state of the 
non-local modes as 



\M,N;0 = 



(Air (Bi 



N 



-m- 



(6) 



B. Two- mode-squeezed number states 

Let us write the annihilation operators of mode A and 
mode B by a = (xa + ipa)/V2 and b = (xb + ipB)/\2. 
This implies the bosonic commutation relations [a, a'] = 
1, [b,^] = 1, and [a, b] = [a, 6 T ] = 0. The number states 
of the local modes are defined by \ti)a '■= (o^)"|0) a/ 'Vn\ 
and \n)B '■= (6^)"|0) b I 'Vn*. together with the conditions 
for the vacuums, a \0)a = and b\0) B = 0. By using 
this standard notation, the two-mode-squeezed (TMS) 



From the construction, it satisfies the orthonormal rela- 
tion (M',N';£\M,N;£) = 5 NiN ,5 m ,m>- It is entangled 
whenever £ ^ and its wave function is non-Gaussian 
except for M = N = 0. From Eqs. ([5]) and © wc 
can confirm the relation \M,N;g) = V^M) A \N) B ■ This 
implies that the TMS number states can be generated by 
acting the two-mode squeezing operation on the product 
of number states as in FIG. QJa). 

The quadrature moments for the TMS number states 
are routinely calculated by using the following relations 
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on the annihilation and creation operators: 



N 
10 r 



— , a 1 



(7) 



The elements of the covariance matrix are determined to 
be 



.4 



D 



C 



l + 2M + (l + 2iV)|£| 2 



\ + 2N +(\ + 2M)\tf 

2\£,\{\ + M + N) r cosO sin0 
1 - |£|2 ^ sin6» - cos6< 



1 

1 

1 
1 



(8) 



where we write £ = |£|(cos# + isin#). This implies an 
explicit formula of the quantity D in Eq. as 



D 



((1 + N)(l + M)\£\ 2 -NM) 



x(MN\£\ 2 -(l + N)(l + M)). 



(9) 



From the condition |£| < 1 we can see that the final 
term is negative, i.e., (MN\£\ 2 - (1 +iV)(l + M)) < 0. 
Hence, the separable inequality D > in Eq. ([3} is 
violated when 



M 



i - iei s 



AT 



< 



iei s 



(i-iei 2 ) 2 ' 



(10) 



The parameter regime of this inequality can be depicted 
in FIG. [2] From Eq. (|10|) . we can specify the subset of 
the TMS number states whose inseparability cannot be 
confirmed by the measurement of the covariance matrix. 
We can see that the majority of the TMS number states 
belong to the class of the purely non-Gaussian entangled 
states. 



C. Beam-splitted number states 

Let us consider the beam-splitter transformation de- 
scribed by 



U r aUl 



Q - rb -., d r :=UM r * a + b 



(11) 



e ratfc r afet is the unitary operator of the 



where U r 

beamsplitter and we assume < \r\ < oo. The field op- 
erators fulfill the bosonic commutation relations [c r , c\] = 
[dr, d\] = 1 and [c r , d r ] = [c r , d\\ = 0. In the case of r = 1 
we have the half-beamsplitter transformation 



a — b a + b 

c = — — , a 



V2 



V2 



(12) 



FIG. 2: A pair of non-negative integers (M, N) assigns a 
TMS number state \M, N; £). In the gray regime, Eq. JTU]) is 
fulfilled, and entanglement of the corresponding TMS num- 
ber states can be verified by the measurement of the covari- 
ance matrix. On the outside (and the boundary) of this 
regime, the covariance matrix of |M, N; £) has to be compat- 
ible with a separable Gaussian state, and its inseparability 
cannot be confirmed by the measurement of the covariance 
matrix. Here, we set the parameter £ = 0.7 in Eq. (|10[) . 



In the limit of \r\ — > 0, there is no interaction between the 
modes and we have a trivial transformation (c, d) = (a, b) 
whereas we have another trivial transformation (c, d) = 
(b, a) in the limit of \r\ — > oo up to the phase factor. 

The beam-splitter transformation produces another 
family of orthonormal states [see, FIG.Q](b)], 



|n, to; r) :=U r \n,m) 



(4) n (4) r 



|0,0). (13) 



We call this state a beam-splitted (BS) number state. 
In the case of n = m — 0, the state |0, 0;r) = |0,0) 
is a separable Gaussian state. Except for this case the 
BS number states are non-Gaussian and entangled as we 
prove in the following. 

Let us assume the non-zero-photon case of n + m > 1 . 
Recall that any Gaussian states are non-orthogonal with 
each other. Since, a non-zero-photon state \n, to; r) is or- 
thogonal to the vacuum state from the orthonormal rela- 
tion (n', to'; r\n, to; r) = 6 n _ n i6 m _ rn i it cannot be a Gaus- 
sian state. From Eqs. (fTTj) and (fTB"]) . we can verify that 
the BS number state |n, to; r) has Schmidt decomposition 
with the Schmidt basis {\n + m — k, fe)}fc=o,i,2,— n+m and 
the Schmidt rank is at most n + to + 1. It is direct to 
see that the coefficients of the edges of the Schmidt ba- 
sis (n + to, 0|n, to; r) and (0, n + m\n, to; r) is non-zero. 
This implies that one can find entanglement from the 
subspace spanned by |n + m, 0) and \0,n + m). Hence, 
the BS number states are non-Gaussian entangled states 
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except for the case of n = m = 0. The inseparability of 
the BS number states can be also proven from the theo- 
rem for beam-splitter entangler [33] (see also [HI)- Note 
that a BS number state is described by a finite num- 
ber of Schmidt bases, and it is thought to be an entity 
of discrete-variable systems. On the contrary, Schmidt 
rank of the TMS number states is infinite, and they are 
thought to be entities of purely infinite-dimensional sys- 
tems. 

The quadrature moments for the BS number states are 
readily calculated by using the following relations on the 
annihilation and creation operators, 

c r + rd r t 4 + r *4 

n — n > — 



^ _ d r - r*c r _ 4 - rc t 



(14) 



The elements of the covariance matrix 7 of Eq. ((2|) are 
determined to be 



A = 
B = 
C = 



1 + \r\ 2 +2(n+ \r\ 2 m) ( 1 



l + |rf V° 1 / ' 

l + |r| 2 + 2(\r\ 2 n + m) ( 1 

l + |rf ^0 1 

2r(m - n) f cos0 -sin</> 



1 -f r 2 V sm V cos ' 



(15) 



where r = \r\(coscj> + isin(j>). From these relations and 
Eq. ©, we have 



D 



1G 



r (m(l + n) + (1 + m)n\r\ 2 ) 



(l + |r| 2 ) 2 ' 
x((l + m)n + m(l + n)|r| 2 ). 



(16) 



This implies that the separable condition D > of Eq. 
§3§ is always satisfied for the BS number states. (This 

result is also given by the fact det C = 4 ^"~|"| 2 |^ > 

with Lemma of Ref. (2^1, that is, Gaussian states with 
det C > are separable.). Hence, they arc classified to 
the purely non-Gaussian entangled states. 

It might be interesting that a relation similar to Eq. 
([TO)) can be derived for the BS number states when we 
consider a higher order entanglement condition given by 
Hillcry and Zubairy (33j : 



(aWb) < |(a6 f )| 2 . 
For the BS number states of Eq. we have 



(17) 



l + |r| 2 / 

x {(1- |r| 2 ) 2 + |r| 2 [m(77i-l)+n(7i-l)]}, 



i&t) 



1+ r 



■(m — n). 



(18) 



Then, the entanglement condition of Eq. (fT7]) turns out 
to be 



\r\ 



•|2 



1 + r 



1 + M 



< 



1 + r 



(19) 



This is an inverse proportional condition on the photon 
numbers of entangled modes for the inseparability similar 
to Eq. (Hig> [See also, FIG. [2]. The condition of Eq. $Jfy 
is satisfied only for the cases of (m,n) = {(fc,0), (0,fc)} 
with k = 1, 2, 3, ■ • • . This is an advantage of the mea- 
surements of higher order moments for detection of the 
purely non- Gaussian entanglement although it can de- 
tect the entanglement of a rather small portion of the BS 
number states. 

Now, it has turned out that the measurement of the co- 
variance matrix is not useful to detect the entanglement 
for a significant portion of the photonic non-Gaussian en- 
tangled states. We thus consider higher order criteria in 
the subsequent section. 



III. VERIFICATION OF PHOTONIC 
NON-GAUSSIAN ENTANGLEMENT VIA THE 
MEASUREMENTS OF HIGHER ORDER 
MOMENTS 

In this section we consider the separable conditions 
based on the operators that form SU(2) algebra and 
SU(1,1) algebra [H [H ■ We show that thc mea- 

surements of quadrature moments to the fourth order are 
sufficient to verify the present non-Gaussian entangled 
states except for a small portion of the BS number states. 
We also show that the entanglement of this small portion 
is detectable by the measurements of further higher order 
moments based on the idea of Ref. [33 ■ 



A. SU(2) and SU(1,1) commutation relations and 
separable conditions 



Let us define the set of operators 



J z 



1 



1 



(<Jb + atf), K x = -(a f tf + ab), 



-(atb-atf), K y 



2i 



(N a ~N b ), K k 



ab), 
N b + 1), (20) 



where N a := a^a and Nb := b^b are number operators of 
the local modes. The operators J; with i = x, y, z form 
the SU(2) algebra while th e op erators Ki with i = x,y,z 
form the SU(1,1) algebra |32l. l33l. |35T - [37| . The commu- 
tation relation [K y ,K z ] = iK x implies the uncertainty 
relation (A 2 K y )(A 2 K z ) > \\{K X )\ 2 . From its partial 
transposition we have the separable condition in Eq. (10) 
of Ref. M 



(A 2 (J,)> + 7 



(A 2 ^) > \\{J X )\\ (21) 
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where the partial transposition of the second sys- 
tem with respect to the number basis is calculated 
by the relation (a tfc a / 6 t °6 p ) PT = (at k a l (tf Wy) = 
(aJ fe a z 6^ p 6°). Similarly, starting from the uncertainty re- 
lation (A 2 J v ) (A 2 J z ) > \\{J X )\ 2 we have another separa- 
ble condition 



(A 2 (^)> - 7 



(A 2 J 2 )>i|(^)| 2 . 



(22) 



In the next two subsections we investigate the role of 
these separable conditions on the detection of the non- 
Gaussian entanglement. The entanglement of the present 
non-Gaussian states is induced by the unitary operators 



V £ = e« a b -«* ab in Eq. © and U r 



(fTT|) , and the Hamiltonians can be connected with a par- 
tial transposition as it maps ra^b—r*ab^ — > ra^b^ — r*ab. 
Hence, it might be insightful to find general relationships 
between the non-Gaussian entanglement and the separa- 
ble conditions with the terms in the Hamiltonian. 



B. Entanglement verification of the TMS number 
states 

The TMS number states of Eq. are simultane- 
ous eigenstates of the number operators andi^. 
From the definition of J z in Eq. ((20)) and the relations 
in Eq. ([5]) we have 



iJif-Bj^ = N a -N b = 2J z 



(23) 



This implies that the TMS number states are eigenstates 
of J z and (A 2 J Z ) = 0. In turn, we can write 



K x = 



1 



2(i-iei 2 r c 

+ (A e + ZBl)(B £ + ZAt)] 



(24) 



Then, its expectation value for the TMS number state 
\M, N; £) of Eq. © is calculated to be 



2(1 -|£| 2 )^ 



(AlAe + BlBt + i) 



2(1 -m 



(M + N + l). 



(25) 



This implies \(K X )\ > 0. Consequently, for the TMS 
number states of Eq. ([6]), the left-hand side of Eq. (|22|) 
is zero and right-hand side of Eq. (|2"2"]l is positive. Hence, 
we can confirm the inseparability of the TMS number 
states from the violation of Eq. (|2"2l . 

Since the separable condition of Eq. (|2"2"]l is violated by 
the TMS number state \M,N;£) for any integers M, N, 
and |£| > 0, it might be valuable to emphasize the mech- 
anism of the entanglement detection. On one hand, the 
TMS number states exhibit strong classical correlation on 
the photon number as we have (A 2 J z ) = 0. On the other 



hand, they show coherence between the (Schmidt) bases 
{\k, k) , \k + 1, k + 1)} as the pair creation and annihila- 
tion of the photons occur coherently. This keeps the term 
(K x ) non-zero. Therefore, we can say that the strong 
classical correlation and the coherence on the Schmidt 
bases signify the existence of entanglement. 

Note that there is a different entanglement criterion 
that can verify the entanglement of the TMS number 
states based on the measurement of quadrature moments 
to the fourth order (2Cj. This criterion is formulated to 
verify the entanglement of the TMS number state when 
the parameter £ is specified. Note also that test of a se- 
quence of high order separable conditions for TMS num- 
ber states have been reported in Ref. [l|| . 



Entanglement verification of the BS number 
states 



The BS number states of Eq. (| 13|) are simultaneous 
eigenstates of the number operators c^c and Sd. From 
the relations in Eq. (fTTI) and the definition of K z in Eq. 
([50"]) we have 



dU = 



a' a ■ 



6 f & = N a + N h = K z - 1. (26) 



Hence, the BS number states are eigenstates of K z , and 
we have (A 2 K Z ) = 0. In turn, from Eq. (fTS]) and defi- 
nition of J x we have (J x ) = ^K^TFF) ( n — m ) ^ or ^ e ^ 
number state \n, m; r) of Eq. (fl3|) . This implies that the 
left-hand side of Eq. ([21]) is zero and the right-hand side 
of Eq. (|2ip is positive whenever n ^ m. Hence, we can 
confirm the inseparability of the BS number state from 
the violation of Eq. (|21[) except for the case of n = m. 

To cover the case of n = m, we consider another sepa- 
rable condition proposed in Ref. [37| . Following Eq. (11) 
of [37j, let us define 



H x — 


i[(ot&t)» 4 


-{abf] 


Hy = 




~ («&) 2 


N + = 


\(N a + N b 


), 


L x = 




(a&t) 2 ] 


Ly = 




(a6t) 2 



(27) 

From the commutation relation [H y ,N + ] = iH x we have 
the uncertainty relation (A 2 H y )(A 2 N + ) > \\{H X )\ 2 . 
Partial transposition on this inequality leads to the sep- 
arable condition in Eq. (14) of [37j . 



<A 2 (4)) + (iV 22 )] <A 2 iV + > > -\{L X )\ 



(28) 



where N 22 := \[a 2 , (a*) 2 ] <g> [b 2 , (6+) 2 ] = (2at a + l)(2&t&- 
1) > 1 is a positive operator. 
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For the BS number states of Eq. (fT3|) we have 
(A 2 N + ) = (A 2 K Z ) = 0, and the left-hand side of Eq. 
([28)) is zero. In turn, by picking the terms which include 
the same number of annihilation and creation operators 
from the expression 

= 2(1+ 1 |r|a)a {[(ct) 2 + (r*rft)='+2r* C t^] 
x[d 2 + {r*c) 2 -2r*cd] 

+ [(dt) 2 + {rJ f - 2rc t d t ][(c 2 ) + r 2 d 2 + 2rcd]} 

(29) 

we have 

{ix) = 2(l+|r|T ((ct)2c2 + (dt)2rf2 ~ 2cfcdU) 
r 2 + (r*) 2 

= 2(1 + |r| 2 ) 2 [m(m ~ ^ + n(?1 ~ ^ ~ 4nm] ' 

(30) 

- .2 _i_ / * \ 2 I 

This implies |(I<a:)| = /i^_L.iaja n{n + 1) > for the case 
of n = m . Hence, the left-hand side of Eq. (f28|) is zero 
and right-hand side of Eq. ([28| is non-zero for the BS 
number states with n = to ^ 0. Therefore, the violation 
of the condition of Eq. (|28p can verify the entanglement 
of the BS number states with n = m. 

Similar to the situation that the TMS number states 
violates the condition of Eq. (|22p , the mechanism of en- 
tanglement detection of the BS number states can be ex- 
plained by the coexistence of strong classical correlation 
and coherence on the photon-number basis. The classical 
correlation is confirmed from (A 2 K Z ) = or (A 2 N + ) = 
while the non-zero coherence of {J x ) > or (L x ) > 
is observed on the (Schmidt) bases with single-photon 
exchange {|n + m — k, k) , \n + m — k — 1, k + 1)} 

or the (Schmidt) bases with two-photon exchange 
{ \n + m — k, k) , \n + m — k — 2, k + 2)} for n ^ m and 
n = to, respectively. Therefore, we can observe basically 
the same structure in non-Gaussian entanglement for 
the TMS number states and the BS number states. 

It might be worth to mention the existence of the BS 
number states that do not violate the higher order sepa- 
rable condition of Eq. (T25|). The set of such BS number 
states is specified by the condition (A 2 iV + ) = 0, namely, 
it can be specified by the pairs of non-negative integers 
(m,n) that fulfill (n— m) 2 — (n+m) = 0. The distribution 
of the pairs (m, n) for 1 < n, to < 10 6 is shown in FIG. [31 
Hence, the separable conditions of Eq. I|28p is also widely 
useful to verify the entanglement of the BS number states 
but there is a small portion of the states which fulfill the 
condition. Consequently, at the present, we require two 
separable conditions in order to verify whole class of the 
BS number states based on the measurement of quadra- 
ture moments. It is in a sharp contrast to the case of the 
TMS number states, in which the entanglement can be 
verified by the violation of the single separable condition 
ofEq. (22). 
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FIG. 3: Distribution of (m, n) that satisfies (n — m) 2 — (n + 
m) = 0. There are eleven solutions for < m < n < 10 6 : 
(0,1), (1,5), (5, 20), (20, 76), (76, 285), (285, 1065), (1065, 
3976), (3976, 14840), (4840, 55385), (55385, 206701), and 
(206701, 771420). The solutions (m,n) = {(0, 1), (1, 0)} are 
not displayed in this figure. 



IV. SUMMARY AND REMARKS 

We have considered two photonic families of entangled 
states generated from the products of number states by 
using two-mode Gaussian entangling operations and in- 
vestigated their inseparability associated with the entan- 
glement criterion based on the measurements of quadra- 
ture moments. We have specified the parameter regime 
of the states where the entanglement cannot be verified 
by the measurement of the covariance matrix, which cov- 
ers the case of the measurement of quadrature moments 
to the second order. We have also shown a couple of 
continuous- variable entanglement criteria to verify whole 
the non-Gaussian entanglement of the two families. An 
interesting finding is that the measurement of the co- 
variance matrices is useless to detect substantial part of 
the non-Gaussian entanglement. On the other hand, it 
has been shown that the measurements to the fourth or- 
der moments are sufficient for the entanglement verifica- 
tion except for a very small portion of the non-Gaussian 
states. 

Generation of non-Gaussian entangled states from the 
number states and Gaussian operation would be an im- 
portant goal of experiments and systematic analysis of 
such states might be an insightful step to comprehend 
quantum correlation beyond the Gaussian framework. 
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